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Gravel beds can prevent sand-dust emission and weaken sand-dust flux. We used wind-tunnel experiments 
and field observations on artificial gravel beds above the Mogao Grottoes to quantifj^ their impact. 
In the report, we identified a significant correlation between gravel roughness and its ability to trap 
wind-transported sand. The optimal combinations of gravel diameter and coverage were determined. The 
greatest roughness is achieved when small gravel coverage is 75%, medium 40% and large 45%. We found 
that initial wind speed and gravel coverage are the key factors controlling the amount of sand trapped by the 
gravel beds. 



The protective measures employed for controlling sand disasters have focused principally on mechanical, 
vegetation and chemical defense\ The Mogao Grottoes (40°5'17.25"N, 94°40'16.67"E) is a famous world 
cultural heritage site located on the desert margins in western Gansu Province, China. They contain 
45,000 m^ of priceless Buddhist murals, 2,500 painted sculptures, and 492 caves^. However, destructive aeolian 
sand has threatened the murals in the Mogao Grottoes since the Five Dynasties (907-960 C.E.). The grotto murals 
have been consistently weathered over the centuries by aeolian sand through wind erosion, wind accumulation 
and dust suppression. Many protective measures have been adopted (Fig. 1), but they have not achieved their 
stated aims^ Prevailing local harmful winds are principally north westerlies, with westerlies and south westerlies of 
secondary importance. When they blow across the land surface, sand accumulation occurs. As a general rule, 
aeolian sand does not blow easily across wind- eroded ^ol?/ surfaces, especially on the borders between gravel beds 
and desert lands. Here, wind-eroded gobi refers to a land surface covered by only a thin layer of large diameter 
gravel, the lower layers being an admixture of sandy gravel and earthy deposits. For hindering the transference of 
sand from nearby Mount Mingsha to the grottoes,creating artificial gravel beds is an effective way to reduce 
aeolian sand damage of the Mogao Grottoes. We developed artificial gravel beds above the Mogao Grottoes with 
varied gravel diameters and coverage based upon the original topography (cobbles, sand gravel beds, and sandy 
gravel beds). Increased surface roughness and resulted in changes to air flow made these artificial gravel beds 
effective for controlling secondary dust flux and weakening the sediment transport intensity of sand- carrying 
wind. 

The earliest studies of gravel beds were by Redway^ and Osborn^. Gobi dust storms are common in arid 
environments^. At the macro-scale, Jugder et aV studied the spatial and temporal distribution of dust mass 
concerntrations of PMIO and PM2.5 in the Gobi Desert in Mongolia, whilst Zhou and Huang^ studied the 
momentum and sensible heat fluxes of turbulent transfer characteristics on a gobi surface in Dunhuang, 
China. These macro-scale studies revealed the characteristics of sand drifts and air flow fields, but studies of 
the concrete interaction between sand and gobi are increasingly the focus of attention. Nickling and Neuman^ 
investigated the deflation drags associated with various diameters and spacings of glass spheres substituted for 
gravel in wind-tunnel tests. The relations between the drag coefficients of gobi surfaces, gravel diameter and 
coverage were defined by Dong et aU^. Subsequently, Mikami et aU^ used a sand particle counter (SPC) with a 
laser beam transmitter and receiver to research saltation fluxes at different heights above ^o^/ surfaces combined 
with meteorological elements. Zhang et aV^ convincingly interpreted the structure of wind-blown sand on gobi 



SCIENTIFIC REPORT; | 4 : 4341 | DO!: 1 0.1 038/srep04341 



1 




Figure 1 | Study site above the Mogao Grottoes. In the gravel bed above the Mogao Grottoes were a sand- accumulation instrument with 16 
orientations (lower left) and a Japanese HOBO weather station (lower right). Upper right is a wind rose derived from wind measurements above the 
Mogao Grottoes during 1992-2008. Area C had a treated surface (large gravel, 30% coverage) and a non-treated surface (medium gravel, 30% coverage); 
Area F had a treated surface (medium gravel, 30% coverage) and a non-treated surface (medium gravel, 30% coverage). The white lines point to the 
measuring sites. The photographs were provided by G. S. Li and W. F. Wang. 



surfaces and found that blown -sand activities on such surfaces are 
mainly concentrated at heights below 20 cm. Gao et a\}^ studied the 
aerodynamic roughness and drag coefficients of gohU whereas Wu^ 
found that gravel beds could increase aerodynamic roughness and 
inhibit wind erosion. 

Generally, the variation of sand weights at different heights 
decreased exponentially^^' ^^ but Qu et aV^ proposed that the distri- 
bution of sand particle weights at different heights on gohi surfaces 
does not simply follow logarithmic or exponential laws. Instead, with 
extreme sand weight the gohi surface height will increase as the wind 
speed increases. Dong et aU^ thought that maximal aerodynamic 
roughness length occurs at gravel coverage ranges between 40%- 
70%. On that basis, Liu et al}^ made a detailed study of the relation 
between gravel coverage and diameter, and concluded that the com- 
bination of 30% of gravel coverage and gravel of a 2-4 cm diameter 



could resist 50% of the sand drift weight. Tan et aV^ proposed that 
the change in the sand transport rates of gohi surfaces in relation to 
height follows a quadratic polynomial. Liu et al and Tan et al in a 
wind-tunnel experiment involving various gravel diameters, gravel 
coverage and sand weights, also demonstrated that by building an 
artificial ^0^/ surface with a gravel coverage of 30%, a dynamic wind- 
blown sand flux equilibrium could be achieved. Zhang et aU^ studied 
the vertical distribution of wind-blown sand flux and attempted to 
quantify sand erosion and deposition above gravel beds. 

At present, few analyses of the effect of optimum gravel diameter, 
gravel roughness and other relative influencing factors on artificial 
sand beds are available. In this study, we have systematically analyzed 
the characteristics of gravel bed sand deposition using field wind- 
tunnel experiments and field observations. In field wind-tunnel 
experiments, the relations between wind speed, abrasion time, sand 
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flux structure, gravel coverage and roughness were quantitatively 
analyzed; and the real impact of sand deposition on gravel beds 
was investigated in the field with sand-accumulation instruments 
and sand- accumulation buckets. 

Results 

Characteristics of the original gravel bed. There is a Japanese 
HOBO weather station located on the gravel bed (Fig. 1) but it is a 
considerable distance from the shifting sands on Mount Mingsha, 
and wind speed and wind direction data (sampling frequency: 
10 min) for 2005-2008 from that station indicate that wind speed 
weakens by the time it reaches the weather station. The Supple- 
mentary Text states that the static pressure threshold wind speed 
of 10.6 m/s at a height of 2 m on the shifting sand surface on 
Mount Mingsha are taken as the threshold wind speed above the 
Mogao Grottoes (see Supplementary Fig. SI online). 

Gravel paving was laid above the grottoes in 2009-2010; even if 
this does not directly affect the threshold wind speed in the shifting 
sand area, the dynamics of this process and the paved area might 
have influenced our analysis, so data from that period was not con- 
sidered in this study. Vis-a-vis the wind direction datasets recorded 
at the HOBO weather station for the purposes of observing gravel 
bed interiors, "sand-blowing wind" is defined as wind with sand 
particles which blows from the N, NNE, NE, ENE, E, ESE, SE, SSE 
or S; "sand-carrying wind," which in fact contains sand particles 
from sand-blowing wind, is defined as the primary harmful wind 
blowing from Mount Mingsha to the grottoes from a SSW, SW, 
WSW, W, WNW, NW or NNW direction. More specifically, sand- 
blowing wind does not carry sand particles when first blowing over 
the gravel bed, but rather collects sand particles as it blows. 
Maximum wind speeds and wind directions for 2005-2008 are 
shown in Supplementary Fig. SI online. Table 1 shows that sand- 
blowing winds and sand-carrying winds most often occur during the 
summer and spring on the original gravel bed above the Mogao 
Grottoes. 

Aeolian sand data collected by sand-accumulation instruments 
with 16 orientations for 2004-2006 (Fig. 1) provide information 
on sand transport characteristics for the original gravel bed above 
the Mogao Grottoes. Supplementary Fig. S2 online presents trans- 
ported sand weight as a percentage of the total from different orien- 
tations as found on the original gravel bed near the grottoes' sides. It 
shows that transported sand weight from a NNE direction was 
20.37% of the total transported sand weight; E, 15.21%; ENE, 
10.73%; NE, 9.63%; ESE, 8.98%; N, 8.14%; SE, 1.97%; S, 1.27%; and 
SSE, 0.93%. The percentages of transported sand weight from a 
NNW direction were 8.21%; WNW, 5.82%; NW, 2.51%; SW, 
1.98%; W, 1.69%; WSW, 1.35%; and SSW 1.20%. Thus, dust storms 
occurring on the original gravel bed are caused principally by sand- 
blowing winds, while the main function of the artificial bed is to 
prevent further sand transport by sand-blowing winds. 

Characteristics of artificial gravel beds. Gravel coverage, wind speed 
and abrasion time. Wind speed and abrasion time determined the 
sand weight deposited on the gravel, and thus ultimately determined 
the coverage of the gravel still exposed above the sand's surface. 



Gravel coverage was also influenced by the roughness of the gravel. 
Also, the sand weight at different heights influenced the coverage of 
the gravel exposed above the sand's surface. 

The relations between wind speed, abrasion time, and coverage of 
small, medium, and large gravels were obtained from 2004 field 
wind-tunnel experimental data. In that experiment, wind speeds 
were quantitatively controlled at about 8 m/s, 10 m/s, 12 m/s, 
14 m/s, 16 m/s and 18 m/s. We used ARC INFO 8.3 software 
(ESRI Corp., Redlands, California) for the vectorization of some 
sample blocks of gravel bed wind-blown for differing times; the 
percentage of total surface area that had gravel >2 mm was then 
entered into the statistical computing functional module, and the 
coverage calculated. The respective relations between wind speed 
Us (m/s) and coverage CoVs (%), and abrasion time (min) and 
coverage CoVs (%) of small gravel can be expressed as: 

=28.169-0.919 X [/, + 0.915 x U^, R^ = 0.99 (1) 

Gov, = 29.770+ 176.838/T,, 7^^ = 0.98 (2) 

The respective relations between wind speed (m/s) and coverage 
CoVm (%), and abrasion time (min) and coverage Cov^ (%) of 
medium gravel can be expressed as: 

Cov^= -18. 755 + 27.644 In [/^, i^^ = 0.98 (3) 

Covm = 73.840 - 3.400T^ + 0. 129T^ - 0.001 T^, = 1.00 (4) 

The respective relations between wind speed Ul (m/s) and coverage 
Covl (%), and abrasion time Ti (min) and coverage Covl (%) of large 
gravel can be expressed as: 

Covi =21.749+ 1.447[7l + 0.093[/^, i^^ = 0.95 (5) 

Covl = 98.904 - 6.218Ti + 0.214T^ - 0.003T^, = 0.94 (6) 

Thus, the relation between wind speed and coverage of small and 
large gravel fits the quadratic function; the relation between abrasion 
time and coverage of medium and large gravel fits the cubic function; 
the relation between abrasion time and coverage of small gravel fits 
the reciprocal function; and the relation between wind speed and 
coverage of medium gravel fits the logarithmic function. The math- 
ematical formats of these factors provide an insight into their 
interactions. 

Sand weight, wind speed and height. Analysis of field wind-tunnel 
experimental data in 201 1 showed that a 20% coverage of small gravel 
approximated to a 23.53% coverage in the original gravel bed. To 
reduce costs, a coverage of 30% of medium and large gravels for the 
artificial gravel bed was shown to be both practical and appropriate. 
As for the underlying sedimentary layer, a coverage of 15%-20% of 
small gravel proved to be adequate. Thus, 20% small gravel and 30% 
medium and large gravels approximated the characteristics of the 
original gravel bed and the artificial gravel bed above the Mogao 
Grottoes, respectively. 

The relation between wind speed Up (m/s) and the correlative sand 
weight Qp (g) at varying heights (1 cm, 3 cm, 5 cm, 9 cm, 13 cm and 
17 cm) under controlled initial wind speeds (8 m/s, 10 m/s, 12 m/s 
and 14 m/s) in the original gravel bed (small gravel, 20% coverage) 
are expressed according to the maximum value of R^ using these 
standard formulae: 

Qp8 = 3.346+ 12.725Up8-3.045[/^8 + 0.180[7^^8, i^^ = 0.96 (7) 
Qpio= -26.338 + 22. 196[/pio-3.448U^^io+0.154[/^\o, i^^=0.99 (8) 
Q^i2= -1.977 + 2. 719Upi2-0.186[7^i2' ^^ = 0.81 (9) 



Table 1 | The frequency and abrasion time for both sand-blowing 
wind and sand-carrying wind (the sand-blowing wind containing 
sand) on the original gravel bed in 2005-2008 

Sand-blowing wind Sand-carrying wind 

Frequency Abrasion time (h) Frequency Abrasion time (h) 

Spring 43 events 28.5 30 events 13 

Summer 54 events 28 51 events 12.7 

Autumn 6 events 26.7 9 events 2.3 
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Qpi4= -1.569+ 1.680[/pi4-0.098[/^\4' ^^ = 0.84 (10) 

where the subscripts p8 through pl4 refer to the initial wind speeds 
described above. These same relations pertaining to a mixture of 
medium and large gravel of 30% coverage are expressed as: 

Qa8 =6.488+ 1.719[/a8-0.460[/^28 + 0.022[/^^8, R^=0.97 (11) 
Qaio = 3.972 + 6.083[7aio - 0.835[7^\o + 0.026[/^\o, = 0.88 (12) 
Qau = - 18.163 + 11.607(7ai2 - 1.322(7^^12 + 0-043(7^^2' =0.87(13) 

Qai4= -47.777+ 19.072[/ai4-1.930[/^^i4 + 0. 0601/^^4, i^^=0.67(l4) 

Thus, initial wind speed (wind speed calculated as net wind when 
blowing into the gravel bed area) and gravel coverage are key factors 
that affect accumulated sand weights in the gravel bed. However, the 
effect of sand weight as related to gravel diameter is not quite as clear- 
cut. The relation between wind speed and accumulated sand weight 
at the same height fits the cubic function at lower initial wind speeds 
in the original gravel bed; this changes into the quadratic function 
with an increase in initial wind speed. However, in the artificial gravel 
bed, regardless of how great the initial wind speed is, the relation 
between wind speed and accumulated sand weight at the same height 
always fits the cubic function, but the correlation becomes smaller 
with an increase in initial wind speed. Understanding more about 
these relations could help us to apply them better in practice. 

According to our analysis of 201 1 field wind-tunnel measurement 
data, sand weight (%) in the original gravel bed (small gravel, 20% 
coverage) presented a certain regularity along a vertical direction 
under the influence of different initial wind speeds. At a height of 
0-6 cm from the surface, the size distribution of sand weight (%) 
were Q (8 m/s)>Q(10 m/s)>Q(12 m/s)>Q(14 m/s); at a height 
of 10-20 cm from the surface, the size distribution of sand weight 
(%) were Q (14 m/s) > Q (12 m/s) > Q (10 m/s) > Q (8 m/s). 
Accumulated sand weight at a height of 0-2 cm was maximal 
(20.95%) when initial wind speed was 8 m/s; sand weight at a height 
of 2-4 cm was maximal (16.82%) at 10 m/s; sand weight at a height 
of 4-6 cm was maximal (14.82%) at 12 m/s; and sand weight at a 
height of 4-6 cm was maximal (14.22%) at 14 m/s (Fig. 2). Thus, the 
sand weight maximum decreased with an increase in initial wind 
speed. 



20-1 
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Figure 2 | Distribution of sand weights at various heights in the original 
gravel bed under different initial wind speeds. Small gravel (20%) bed is 
very close to the coverage 23.53% of the original gravel bed. 



Sand weight outlines in medium gravel (30% and 40%) were more 
regular than in large gravel (30% and 40%). Vis-a-vis the size distri- 
bution of sand weight in medium gravel (30%) at a height of 8- 
20 cm, it was Q (14 m/s) > Q (12 m/s) > Q (10 m/s) > Q (8 m/ 
s) (Fig. 3al); in large gravel (30%) at heights of 8-20 cm and 14- 
20 cm it was Q (14 m/s) > Q (12 m/s) > Q (10 m/s) > Q (8 m/s) 
(Fig. 3bl); in medium gravel (40%) at a height of 6-10 cm it was Q 
(8 m/s) > Q (10 m/s) > Q (12 m/s) > Q (14 m/s) (Fig. 3a2); and in 
large gravel (40%) at a height of 18-20 cm it was Q (14 m/s) > Q 
(12 m/s) > Q (10 m/s) > Q (8 m/s) (Fig. 3b2). Thus, the distribution 
of accumulated sand weights at increasing heights became irregular 
with increases in gravel diameter and coverage. 

This demonstrates that accumulated sand weight in the original 
gravel bed follows a definite law along a vertical direction in response 
to different initial wind speeds. Sand weight decreased with an 
increase in initial wind speeds at a height of 0-8 cm, and increased 
with an increase in initial wind speeds at a height of 10-20 cm. The 
distribution of accumulated sand weights in medium gravel at vari- 
ous heights was more regular than that in large gravel on the artificial 
gravel bed. 

Roughness and gravel coverage. Our analysis of 2011 field wind-tun- 
nel experimental data indicated that roughness associated with dif- 
ferent gravel coverage can be revealed by an equilibrium wind speed 
profiler when the height of the wind speed profiler and the net wind 
speed are combined with Supplementary Text Eq. (16). The rough- 
ness curves of small, medium and large gravels presented four peaks 
related to the increase in coverage (Fig. 4). Fig. 4 suggests that the 
following combinations of diameter and coverage are appropriate for 
gravel paving on the gravel bed above the grottoes: small gravel, 75%; 
medium gravel, 40%; and large gravel, 45%. To reduce the labor costs 
of gravel paving the area, the small gravel could be the underlayment, 
with an emphasis for overlay on the medium and large gravels. The 
protective effect of gravel paving can thus be maximized by following 
this construction sequence: first, small gravel (1-2 cm) is spread to a 
ground coverage up to 15%-20% on a treated earth surface, and then 
larger gravel (3-5 cm) is randomly spread to a ground coverage up to 
30%-40%. 

Comparison of original and artificial gravel beds. Field wind- 
tunnel experiments. The sand accumulation rate is expressed as 
sand weight deposited from sand-carrying wind per unit width 
against unit time. In our study, the size distribution of sand 
accumulation rates on the original gravel bed (small gravel, 20% 
coverage) at a height of 0-14 cm was Qd (14 m/s) > Qd (12 m/s) 

> Qd (10 m/s) > Qd (8 m/s) > Qd (6 m/s) (Fig. 5). The sand 
accumulation rate thus increased with an increase in initial wind 
speed: the higher the initial wind speed blowing into the original 
gravel bed area, the greater the rate of sand accumulation. 

The impact of sand accumulation in medium (30% coverage) 
(Fig. 6al) and large gravel (30%) (Fig. 6b 1) areas in the artificial 
gravel bed was greater than in the original gravel bed (small gravel, 
20%) (Fig. 5). The size distribution of sand accumulation rates in the 
medium gravel bed with 30% coverage at a height of 0-14 cm was Qd 
(14 m/s) > Qd (12 m/s) > Qd (10 m/s) > Qd (8 m/s) > Qd (6 m/s). 
Sand accumulation rates at height 0-2 cm were maximal at initial 
wind speeds of 6 m/s, 8 m/s, 10 m/s and 12 m/s, but the sand accu- 
mulation rate at height 6-8 cm was maximal with an initial wind 
speed of 14 m/s (4.65 g/(min-m)) (Fig. 6al). The sand accumulation 
rate over large gravel (30% coverage) at height 0-2 cm was up to 
7.34 g/(min*m) with an initial wind speed of 12 m/s (Fig. 6bl). The 
size distribution of the sand accumulation rate of medium gravel 
(40% coverage) at height 0-12 cm was Qd (14 m/s) > Qd (12 m/s) 

> Qd (10 m/s) > Qd (8 m/s) > Qd (6 m/s) (Fig. 6a2). The size 
distribution of sand accumulation rates for large gravel (40% cov- 
erage) were more irregular (Fig. 6b2). 
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Figure 3 | Distribution of sand weights at various heights in the artificial gravel bed under different initial wind speeds, (al) Variations in 
accumulated sand weights at various heights on a medium gravel (2-4 cm) bed with coverage of 30% under different initial wind speeds (8 m/s, 10 m/s, 
12 m/s, and 14 m/s). (a2) On a medium gravel (2-4 cm) bed with coverage of 40%. (bl) On a large gravel (4-6 cm) bed with coverage of 30%. (b2) On a 
large gravel (4-6 cm) bed with coverage of 40%. 



On the artificial gravel bed, the sand accumulation rate for small 
gravel and its consequent impact became greater with increased 
coverage. Conversely, the sand accumulation rate for large gravel 
became smaller with increased coverage and its overall impact 
became negligible. Therefore, on the whole, the impact of sand accu- 
mulation for a combination of small (40% coverage) and large gravel 
(30%) proved optimal. In addition to the gravel paving's main func- 
tion of preventing sand-driven by wind, intercepting sand particles 
from sand- carrying wind also serves a vital purpose. 

Field observation experiments. We analyzed a sand accumulation 
section (W E) on the northern side of a V-shaped nylon net^ in 
Areas C (treated surface: large gravel, 30%; non-treated surface: med- 
ium gravel, 30%) and F (treated surface: medium gravel, 30%; non- 
treated surface: medium gravel, 30%) of an experimental area of 
gravel paving (Fig. 1). Data are for the periods March 18 - April 
15, 2009 (Phase 1) and May 13 - May 17, 2010 (Phase 2). The 
purpose was to study the effect of sand accumulation on original 
and artificial gravel beds (large gravel 30%, medium gravel 30%). 

In Phase 1, the easterly wind had an initial wind speed of more 
than 10.6 m/s at a height of 2 m (average wind speed was 1 1.67 m/s, 
and abrasion time 30 min); this was defined as the net wind speed 



(sand-blowing wind). In Phase 2 the westerly wind had a wind speed 
of more than 10.6 m/s at a height of 2 m (average wind speed was 
11.77 m/s, abrasion time 15 min), so this was defined as the sand- 
carrying wind. In Phase 1, the average sand accumulation rate over 
the original gravel bed was 0.042 g/(min*m) when a sand-blowing 
wind was present; in Phase 2, the average was —0.171 g/(min*m) 
after the sand-carrying wind had passed. In Phase 1, average sand 
accumulation rates for large (30%, treated surface) and medium 
gravel (30%, treated surface) were —0.219 g/(min-m) and 6.075 g/ 
(min*m), respectively, when a sand-blowing wind was present; In 
Phase 2, these averages were 0.292 g/(min-m) and 1.938 g/(min-m), 
respectively, after the sand- carrying wind had passed. 

Contrary to expectations, over the original gravel bed the sand 
accumulation rate was positive when a sand-blowing wind was pre- 
sent, and negative after a sand-carrying wind had passed. These 
abnormalities in sand accumulation rates in the original gravel bed 
are an important cause of severe aeolian sand damage in the Mogao 
Grottoes. The sand accumulation rate of large gravel (30% coverage, 
treated surface) in the paved gravel experimental area was negative 
when a sand-blowing wind was present and positive after a sand- 
carrying wind had passed; this shows that large gravel (30%, treated 
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Figure 4 | Variations in the roughness of small, medium, and large gravel 
against increases in coverage. Coverage of relatively small gravel (and 
aerodynamic roughness lengths) are 25% (0.185 cm), 45% (0.354 cm), 
55% (0.352 cm) and 75% (0.471 cm); those of medium gravel are 20% 
(0.193 cm), 40% (0.366 cm), 50% (0.300 cm) and 60% (0.303 cm); and 
those of large gravel are 20% (0.211 cm), 30% (0.379 cm), 45% 
(0.407 cm) and 75% (0.265 cm). 

surface) corrected the abnormal sand accumulation rate in the ori- 
ginal gravel bed. The sand accumulation rate of medium gravel (30%, 
treated surface) in the paved gravel experimental area was positive 
when a sand-blowing wind was present, similar to that of the original 
gravel bed, but the sand accumulation rate was also positive at 
1.938 g/(min-m) after a sand-carrying wind had passed; this rate 
was significantly greater than the sand accumulation rate of 
0.292 g/(min*m) for large gravel (30%, treated surface). This sug- 
gests that the impact of sand accumulation on medium gravel (30%, 
treated surface) is more effective vis-a-vis protection against sand 
than on large gravel (30%, treated surface) when a sand- carrying 
wind passes. Therefore, the data suggest that gravel paving using a 
combination of large gravel (30%) and medium gravel (30%) is both 
appropriate and effective. 

On June 2, 2011 a great dust storm (NW wind) occurred in the 
vicinity of the grottoes, and measurement data were acquired from 
sand- accumulation instruments stationed at a height of 1 m along a 
NW SE orientation above the grottoes. The weight of accumulated 
sand in the paved gravel area was approximately 21% less than in the 
shifting sand area at height 0-18 cm (Fig. 7). This indicates that the 
gravel paving captured 21% of the sand from the sand-carrying wind, 
with effective sand accumulation rates. This gravel paving model 
thus shows significant apphcability. 

Discussion 

The gravel paving of the original gravel bed above the Mogao 
Grottoes addresses many research questions relating to gravel 
armoring to reduce wind erosion. Its protective mechanism has sci- 
entific validity and plays an important role in aeolian sand protection 
for infrastructure such as railways, spaceflight, satellite and military 
bases in aeolian sand and gobi areas. 

In the relation between gravel coverage and gravel roughness, we 
determined optimal coverage by gravel with different diameters 
using the peak values of gravel roughness. This fundamentally agreed 
with the perspective of Dong et aV^ that aerodynamic roughness 
length possesses maximal value when gravel coverage is in the range 
of 40%-70%. However, we further refined the specific optimal gravel 
coverage, diameter and roughness, as shown in Fig. 4. The peaks in 
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Figure 5 | Distribution of sand-accumulation rates at various heights on 
the original gravel bed under different initial wind speeds. Sand- 
accumulation rates (g/(min-m)) of a small gravel (20%) bed under 
different initial wind speeds (6 m/s, 8 m/s, 10 m/s, 12 m/s and 14 m/s). 

Fig. 4 indicate that small gravel has the greatest roughness at —75% 
coverage; medium gravel has the greatest roughness at —40% cov- 
erage; and large gravel has the greatest roughness at —45% coverage. 
We proposed that sand particles in the gravel bed would be blown by 
a sand-blowing wind and deposited by a sand-carrying wind. This 
differentiation between the blowing of unsaturated sand drift and the 
accumulation of saturated sand drift in the gravel bed is in accord- 
ance with the findings of Zhang et aP°. It also showed that the sand 
accumulation rate could be regarded as an index for efficiently mea- 
suring erosion and deposition in the gravel bed. 

However, the secondary flow and size effect problems caused by 
the wind tunnel wall in our study remained unresolved. In observa- 
tions of aeolian sand transport, a field wind tunnel could ideally be 
used to simulate the pulsation of wind speed and wind direction 
resulting from thermal convection caused by actual ground surface 
temperature changes. Such a tunnel could also be used to model dust- 
rolling phenomena. Neither of these were possible with the wind 
tunnel that was available to us. In future, we intend to use a com- 
bination of a particle image velocimetry (PIV) system and a field 
wind tunnel in the gravel bed to determine how gravel characteristics 
influence the incidence of sand particles and the distribution of fluc- 
tuating speeds. Problems were also posed by density and 3-D velocity 
distributions for sand saltation rates vis-a-vis different gravel dia- 
meters, coverage, and shapes at different heights in the gravel bed. 

Other drawbacks in our study included the automatic weather 
station not being set up on the surface of shifting sand to monitor 
change when measuring static pressure threshold wind speed for 
shifting sand surfaces. There were pulses in the wind speed from 
the field wind-tunnel fan, indicating that it was not controlled accur- 
ately; errors were in the range of 0.3 m/s. Our calculation of the 
weight of sand carried by the sand-carrying wind as reduced by 
gravel paving was based on only one dust storm where it should have 
been calculated from the average of many dust storms. Further, 
human and natural factors were not accurately controlled, and thus 
possibly disturbed the data. 

Methods 

Gravel bed processing. After proposing the artificial gravel bed concept we 
performed some experiments. Our experimental method was quite complicated, 
involving three gravel diameters and different coverage (15%, 30%, 30%-40%, 45%, 
and 50%-60%). First, the study surface had to be treated by spraying water on the 
earthen surface to create a wet sand layer up to about 5 cm thick, then a stone roller 
was used to grind the surface to make it level. Broken clay was evenly spread on the 
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Figure 6 | Distribution of sand accumulation rates at various heights on the artificial gravel bed under different initial wind speeds, (al) Variations in 
sand accumulation rates (g/(min*m)) at various heights on a medium gravel (2-4 cm) bed with coverage of 30% under different initial wind speeds 
(6 m/s, 8 m/s, 10 m/s, 12 m/s, and 14 m/s). (a2) On a medium gravel (2-4 cm) bed with coverage of 40%. (bl) On a large gravel (4-6 cm) bed with 
coverage of 30%. (b2) On a large gravel (4-6 cm) bed with coverage of 40%. 



surface (each clay batch weighed about 0.5 kg/m^) and a broom was used to mix the 
surface layer clay and the sand soil. The stone roller ground the surface level again. 
Then small gravel (diameter 0.5-2 cm) was evenly spread on the leveled earth surface, 
with a coverage of about 15%. The stone roller was used to compact the admixture 
again to increase earth surface hardness to an average measured value of > 12 N/mm^ 
as verified by a soil hardness tester. Then small gravel (diameter 1-2 cm) was once 
again evenly spread on the prepared earth surface, with a coverage of about 15%. 
Finally, gravels of different diameters were evenly spread on the treated earth surface 
at different coverage. For fair testing and equivalence, gravels of differing diameters 
were directly and randomly tossed onto non-treated surfaces. 

Gravel paving to reduce aeolian sand transport and damage above the Mogao 
Grottoes began in 2009-2010. Our main study site was the typical original gravel bed 
topography and the artificial gravel bed formed after the installation of that paving. 
The original gravel bed topography was made up of cobbles, a mixed sand/gravel bed, 
and a sandy gravel bed. Gravel paving and sand covered 686,650 m^, with coverage 
after ground treatment measured at 991,000 m^, giving a total area of 1,677,650 m^, 
approximately equivalent to the contiguous area of the original gravel bed. An area of 
shifting sands on nearby Mount Mingsha provided another small study block. 

Wind tunnel experiments. Our gravel diameter sampling areas were set at 900 cm^. 
We sampled from the four corners, the midpoints of the four sides, and the cross- 
point of two diagonals, producing a total of nine sampling spots. In each 900 cm^ 
sampling plot, the diameters of 100 randomly selected gravels were measured and the 
average coverage in the original gravel bed calculated. Supplementary Fig. S3 online 
shows that the paved gravel on the original gravel bed was predominantly of small 
gravel (<2 cm). In 2004 the northern experimental area above the Mogao Grottoes 
had gravel coverage of 23.53% in its original gravel bed. Changes in sand coverage on 



gravel of different diameters, at different wind speeds and for different abrasion time, 
was studied based on field wind tunnel experimental data from 2004. Functional 
relations between gravel coverage, wind speed, and abrasion time were defined in 
2004 field wind-tunnel experiments using small (average diameter 0.89 cm), medium 
(average diameter 2.68 cm) and large gravel (average diameter 5.49 cm). 

The threshold wind speed (defined below) of aeolian sand above the grottoes was 
validated by field wind-tunnel (See Supplementary Fig. S4 online) experimental data 
extracted from the surface of a shifting sand area on Mount Mingsha in 2010. Wind 
speed outlines, sand transport, and sand accumulations in clean wind and sand- 
carrying wind conditions in gravel paving blocks with different gravel diameters and 
coverage were obtained from field wind-tunnel data collected in 2011. 

Box-like sand-accumulation instruments (See Supplementary Fig. S5 online) of 
height 20 cm used in field wind tunnel experiments were set up at 3.4 m intervals on 
the gravel bed, spanning different gravel coverage and diameters. Weighed-out sand 
quantities (g) at different heights and the sand accumulation rates (g/(min-m)) on the 
gravel bed were calculated for different lengths of sand-carrying wind blowing times. 
A wind speed profiler with a pitot tube (See Supplementary Fig. S6 online) was also 
used. 

Field experimental observations. Data pertaining to historical sand-accumulation 
conditions were obtained from sand-accumulation instruments with 16 orientations 
(Fig. 1) during 2004-2006. These reflected sand transport characteristics on the 
original gravel bed that had not been destroyed by paving. The transported sand, once 
collected, was weighed by electronic scales with +/— 0.001 g precision. 

During 2005-2008, data pertaining to wind speed and wind direction in the ori- 
ginal gravel bed were obtained from the Japanese HOBO weather station (Fig. 1). By 
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Figure 7 | Distribution of sand weights at various heights in the 
shifting sand area and the gravel paving area above the Mogao Grottoes 
during the great dust storm of June 2, 2011. The wind direction of the 
storm was northwest. Staged sand accumulation instruments of height 
1 m were set up on the shifting sand area and the artificial gravel area along 
a NW SE orientation. 

combining field wind-tunnel experimental data and mathematical formulae, we 
extracted both wind speeds and wind directions. 

This study reports the results from Area C in Fig. 1 (treated surface: large gravel, 
30%; non-treated surface: medium gravel, 30%) and Area F (treated surface: medium 
gravel, 30%; non-treated surface: medium gravel, 30%). Wind speed and wind dir- 
ection were extracted from automatic weather station data correlated with measured 
time intervals for sand weight in the sand-accumulation buckets. The local gravel bed 
sand-blowing and sand-carrying winds were assessed using data from the automatic 
weather station on above the grottoes for the dates March 18 - April 15, 2009 and May 
13 - May 17, 2010. Sand weights collected by sand- accumulation buckets^^ in paved 
gravel experimental blocks with different treatments (dependent upon coverage, 
diameter, and construction technology) during these periods were used to study the 
impact and behavior, under natural conditions, of sand accumulation in paved gravel 
areas possessing various combinations of different gravel diameters and coverage. 1 1 
sand-accumulation buckets were spaced 100 m apart on the northern side of a V- 
shaped nylon net (there were sand accumulation buckets in the original gravel bed 
and because buckets 10 and 1 1 were very close to the net and vulnerable to its effects, 
they were weeded out)^*^. The sand was weighed by electronic scales with +/ — 1 g 
precision. 

The vertical heights of aeolian sand saltation and sand transport under natural 
conditions were measured using sand-accumulation instruments of height 1 m. 
These were set up on the surface of the shifting sand and artificial gravel bed areas 
along a NW SE orientation above the grottoes during a large dust storm (northwest 
wind) on June 2, 2011. These staged sand- accumulation instruments (with openings 
at 2 cm increments) measured the intensity of the dust storm; electronic scales with 
+ /— 0.01 g precision were used to measure sand weights from different heights. 
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